Preprint typeset in JHEP style - HYPER VERSION 



CERN-TH/2001-225 
DAMTP-2001-51 



Brane-World Cosmology, Bulk Scalars and 
Perturbations 



Philippe Brax* 

Theoretical Physics Division, CERN 
CH-1211 Geneva 23. 
philippe . braxOcern . ch 

Carsten van de Bruck & Anne C. Davis 

DAMTP, Centre for Mathematical Sciences 

Cambridge University, Wilberforce Road, Cambridge, CBS OWA, UK. 
C . VanDeBruckOdamtp . cam. ac.uk, a. c . davisOdamtp . cam. ac .uk 



Abstract: We investigate aspects of cosmology in brane world theories with a bulk scalar 
field. We concentrate on a recent model motivated from supergravity in singular spaces. 
After discussing the background evolution of such a brane-world, we present the evolution 
of the density contrast. We compare our results to those obtained in the (second) Randall- 
Sundrum scenario and usual 4D scalar-tensor theories. 



Keywords: supergravity models, cosmology of theories beyond the SM, physics of the 



early universe, large scale structure formations 



*On leave of absence from Service de Physique Theorique, CEA-Saclay, F-91191 Gif/Yvette cedex France 



Contents 



i 
I 



I 



I 



Introduction 

Brane Cosmology 
|2.l| The Background Evolution 
Radiation Dominated Eras 
Matter Dominated Era 
Supergravity Era 

The Matter-Supergravity Transition 



2.2 



Cosmological Perturbations using the Fluid Flow Approach 
|3.1| The Raychaudhuri and Klein-Gordon Equation 



3.2| The Evolution Equation for E^v 
3.3 The Perturbed Dynamics 



Time Evolution of Cosmological Perturbations 

!T| The Randall-Sundrum scenario 
|4.1.l| The high energy era 
|4Ll The low -energy radiation era 



4.2 



4. 1 .3| The matter dominated epoch 

Effects of the bulk scalar field on cosmological perturbations 



4.2.l| The high-energy regime 

4.2.2 Radiation domination 

4.2.3 Matter domination 



4.2.4 Observational consequences 



|fj|. Conclusions 



Supersymmetric Backgrounds 
[A.l BPS configurations 



A. 2 Breaking Conformal Flatness 







10 

n 

u 

16 

n 

M 

18 

1 
1 

1 

19 

m 



23 

24 
H 

m 



1. Introduction 

The idea that we live on a hypersurface embedded in a higher dimensional space has sparked 
a lot of interest recently. These theories are motivated from string theory, where higher- 
dimensional objects, such as D-branes, play an essential role [|lj. Similarily, compactifying 
M-theory (or its effective low-energy limit: 11-D supergravity) on a Si /^-orbifold and 
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compactifying six dimensions on a Calabi-Yau manifold, results in a five-dimensional brane 
world scenario with two hypersurfaces, each located at the orbifold fixed points (see e.g. 
§ and §). 

Brane world theories predict that our universe was higher-dimensional in the past. 
Because of this, there is the hope that certain questions which cannot be answered within 
the context of the standard model of cosmology, can be addressed within these theories. 
Furthermore, cosmology should be a way to severely constrain parameters in these models. 

So far, most cosmological considerations of brane worlds centered around the one- 
brane scenario of Randall and Sundrum Q. In this model, the three-dimensional brane 
universe is embedded in a 5-dimensional Anti-de Sitter (AdS) spacetime. In particular, 
the bulk-space is empty, the only contribution to the curvature comes from the negative 
cosmological constant in the bulk. This simple model already leads to new effects which 
are interesting for cosmology [||. However, most scenarios motivated from particle physics 
predict matter in the bulk, such as scalar fields. In 5D heterotic M-theory, for example, 
one particular scalar field measures the deformation of the Calabi-Yau manifold, on which 
six other small dimensions are compactified |3]]. Other models, motivated by supergravity 
(SUGRA), also predict bulk matter, whose form is dictated by the field theory under 
consideration. 

Cosmology may be a fruitful field where the above ideas can be tested. As such, the 
study of the evolution of cosmological perturbations is extremely important ||, because 
the higher-dimensional nature of the world can leave traces in the distribution of matter 
and/or anisotropies in the microwave sky. A lot of papers investigated different aspects of 
perturbations in brane world scenarios 

The aim of this paper is to develop an understanding of the evolution of perturbations 
in brane world scenarios, in which scalar field(s) are present in the bulk. As a toy model, 
we will use a cosmological realization of a supergravity model in singular spaces M. The 
evolution of the brane world was discussed in depth in || and [[!(]] (see also [jll]] for a 
discussion on brane cosmology and bulk scalar fields). In particular four different cosmo- 
logical eras have been identified in this model. At high energy above the brane tension the 
cosmology is non-conventional before entering the radiation epoch where the scalar field 
is frozen. After matter-radiation equality the scalar field starts evolving in time leading 
to a slow-down of the expansion rate compared to FRW cosmology in the matter domi- 
nated era. Eventually the scalar field dynamics becomes the dominant one leading to a 
supergravity era. Requiring that coincidence between the matter and scalar field energy 
densities occurs in the recent past leads a cosmological constant with a fine-tuning of the 
supersymmetry breaking tension on the brane. In the supergravity era it has been shown 
that the observed acceleration of the expansion of the universe can be understood, i.e. the 
computed acceleration paremeter qo = —4/7 is within the experimental ball-park. This 
accelerated expansion is driven by a bulk scalar field, whose parameters are constrained by 
the gauged supergravity theory in the bulk. On top of this, the model predicts a significant 
evolution of the induced gravitational constant G on the brane. Indeed, the value of G can 
be seen to have changed by 37 percent since radiation/matter equality. Thus, the study of 
cosmological perturbation theory in these models is rather important, in particular the time 
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evolution of G may leave traces in the evolution of perturbations. In this paper we discuss 
the evolution of the density constrast and the effect of the bulk scalar field on cosmological 
perturbations. In Section 2 we review Einstein's equations induced on our brane world and 
discuss the Friedmann equation. We also review the background solutions found in [10], 
which are needed in order to solve the perturbation equations. In Section 3 we derive the 
perturbation equations using the fluid flow approach. This approach is very transparent 
for our purposes and makes it easy to derive the necessary evolution equations. In Section 
4 we discuss some solutions of these equations and discuss their properties. We point out 
the differences to the Randall-Sundrum scenario and usual four-dimensional scalar-tensor 
theories. We conclude in Section 5. In the appendix we discuss some details concerning 
the issue of supersymmetry breaking and conformal flatness. 



2. Brane Cosmology 

In this section we discuss the field equations on the brane and discuss the background 
evolution. 

2.1 The Background Evolution 

We consider our universe to be a boundary of a five dimensional space-time. The embedding 
is chosen such that our brane-world sits at the origin of the fifth dimension. We impose a 
Z2 symmetry along the fifth dimension and identify x§ with — X5. Our brane-world carries 
two types of matter, the standard model fields at sufficiently large energy and ordinary 
matter and radiation at lower energy. We also assume that gravity propagates in the bulk 
where a scalar field eft lives. This scalar field couples to the standard model fields living on 
the brane-world. At low energy when the standard model fields have condensed and the 
electro-weak and hadronic phase transitions have taken place, the coupling of the scalar 
field to the brane-world realizes the mechanism proposed in [|l^] with a self-tuning of the 
brane tension. In this section we derive the brane cosmology equations describing the 
coupling between ordinary matter on the brane and a scalar field in the bulk. 
Consider the bulk action 

Sbulk = / d5x V^(R - lm) 2 + U(<P))) (2.1) 
where n\ = 1/Mg and the boundary action 

Sb = - A / d A xyf^iUB{4>v) (2.2) 
where 0o is the boundary value of the scalar field. The Einstein equations read 

Gab = Rab ~ -^R§ab = T a b + S X5 T^ b (2.3) 

where T ab is the bulk energy-momentum tensor and Tr? is the boundary contribution. The 
bulk term is 

Tab = \ (d a 4db<f> - lg a b(dcf>) 2 ) - ^g ab U (2.4) 
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and the boundary term 



Tg = -\ 9 a b U B ^) (2.5) 

with a, b = ... 3 in the last equation. Following the self-tuning proposal we interpret Ub 
as arising from a direct coupling Ug to the brane degrees of freedom, i.e. the standard 
model fields The vacuum energy generated by the $ l 's yields the effective coupling 

=< V(*) > U% (2.6) 

where the dimension four potential V(&) represents all the contributions due to the fields 
after inclusion of condensations, phase transitions and radiative corrections. 
We also consider that ordinary matter lives on the brane with a diagonal energy mo- 
mentum tensor 

T a b = dmg(-p,p,p,p) (2.7) 

and an equation of state p = up. 

We will be mainly concerned with models derived from supergravity in singular spaces. 
They involve N = 2 supergravity with vector multiplets. When supergravity in the bulk 
couples to the boundary in a supersymmetric way the Lagrangian is entirely specified by 
the superpotential 

U B = W (2.8) 

and the bulk potential 
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If one considers a single vector supermultiplet then supersymmetry imposes that 

W = fe Q * (2.10) 

where a = l/\/3, — l/vl2, these values arising from the parametrisation of the moduli 
space of the vector multiplets, and £ is a characteristic scale related to the brane tension. 

Since supersymmetry is not observed in nature, one should incorporate supersymmetry 
breaking. A natural way to break supersymmetry is by coupling the bulk scalar field to 
brane fields fixed at their vevs. This leads to 

U B = TW (2.11) 

where T = 1 is the supersymmetric case. Larger values of T correspond to supersymmetry 
breaking effects with a positive energy density on the brane. We will analyse the dynamics 
of the coupled system comprising gravity, the scalar field <j> and matter on the brane. 

The dynamics of the brane world is specified by the four dimensional Einstein equations 



Gab = -^-riab + —rT ab + TT ab + ]-V a (fN b4> ~ 777 (V</>) 2 ™a& ~ E a b, (2.12) 

8 4 2 lb 

where r ao is the matter energy momentum tensor, n a {, the induced metric on the brane and 
E ao the projected Weyl tensor of the bulk onto the brane. It appears as an effective energy 
momentum tensor on the brane, called the Weyl fluid ||. 
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The tensor ir a b is quadratic in the matter energy momentum tensor 

T T acT b T a b c d 2\ r n -i o\ 

Kab = -r^T ab ~ + ^y 3TcdT ~ T ) V 2 - 13 ) 

where r = r a a , The Bianchi identity V a G ao = leads to the conservation equation 

V a E ab = -^r ab + V a ir ab + V a P ab , (2.14) 

where the tensor P ab is defined by 

Pab = ~^-n ab + \v a <jH h <t> - ^(V0) 2 n ab (2.15) 
8 2 lb 



This consistency equation ( 2.14 ) will allow us to follow the time evolution of £"00 in the 
background and at the perturbative level. This is crucial in order to define the time 
evolution of the background and of the scalar field and matter perturbations. 

The dynamics of the scalar field is specified by the Klein-Gordon equation which reads 



.0 , rdU B dV 



v * + s-*= (2 ' 16) 

where we have defined the loss parameter by 

and d 2 4>\o stands for the second normal derivative of the scalar field at the brane location. 
The effective scalar potential is defined by 

V= -. (2.18) 

We derive these equations in section 3. 

In order to illuminate the role of Ub further, we point out that the projected Klein- 
Gordon equation can be seen as an equation for the brane energy-momentum tensor 

T ab = V a 0V 6 - ^((V0) 2 + 2V)n ab , (2.19) 

which reads 

V a T ab = ~V b U b - A$ 2 V 6 (2.20) 
o 

and can be derived from the four dimensional effective action for (j> 

S = f d*xy/=gi ((V0) 2 + 2V- T -U B - 2A$ 2 0) , (2.21) 

whenever A<I> 2 is constant. It is remarkable that Ub plays the role of a dilaton coupled to 
the trace of the matter energy- momentum tensor. 

For the background cosmology the Klein-Gordon equation reduces to (see section 3) 

4> + 4H^- T -U' B = -V + A&2 (2.22) 
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where dot stands for the proper time derivative. The bulk expansion rate is defined by 
4H = d r In y/—g\o evaluated on the brane. 

The background cosmology is characterized by its isotropy so we consider that Eoi = 
0, Eij = 0, i ^ j. Moreover we assume that we obtain a FRW induced metric on the brane. 
Due to the tracelessness of E a b it is then sufficient to obtain the differential equation for 

3 ; 2 . 3 \ 3 ; 2 , u b 



Eqo + 4H B E 0Q = dt[—^ + -V\ + -H B tf + -^p, (2.23) 



where Hb is the brane expansion rate 2>Hb = d T In ^—g B \o- This leads to 

E 00 = 1 J dta* ^ (A^ 2 + |v) + j^ 2 + ^ . (2.24) 
Writing the bulk metric as 

ds 2 = e 2A(z,t b ) { _ dt 2 + + e 2B( Z ,t b ) dx 2 t ^ 

the proper time on the brane is defined by dt = e A ^' tb ^dtb- We can always choose the 
boundary condition A(0,tb) = B(0,tb) in such a way 

H = H B . (2.26) 

Together with the Klein-Gordon equation this leads to 

In particular we find that conformal flatness is broken as soon as <ft becomes time-dependent 
and either matter is present on the brane or the energy loss parameter does not vanish. 
Notice that a sufficient condition for breaking conformal flatness is that Newton's constant 
becomes time-dependent. We will comment on this expression when discussing the various 
cosmological eras. 
Using 

Goo = 3-ffg (2.28) 
and after one integration by parts we obtain the Friedmann equation 

rr2 P 2 U B If, da 4 . 1 f , dU B C . 

H l = ^ + ^P + —Jdt—(2V-^)-—J dt^p+-, (2.29) 



36 12 r 16a 4 / dt y r ' 12a 4 / dt 



This equation has already been derived in [|T(]] . Notice that the scalar field <j) enters in the 
definition of the effective Newton's constant. The last three terms are a combination of the 
energy flow onto or away from the brane and the changes of the pressure along the fifth 



dimension [14|. 



Using the conservation of matter 

V a r fc a = (2.30) 
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we find that 

p=-3H B (p + p). (2.31) 

This completes the description of the three equations determining the brane cosmology, 
i.e. the Klein-Gordon equation, the Friedmann equation and the conservation equation. 

Notice that there are two entities which depend on the bulk. First of all there is the 
dark radiation term C/a 4 whose origin springs from the possibility of black-hole formation 
in the bulk. Then there is the loss parameter A<3?2 which depends on the evolution of the 
scalar field in the bulk. It specifies the part of the evolution of cf>, which is not constrained 
by considerations of the brane dynamics. 

In the following we shall describe the case where the bulk theory is N = 2 supergravity 
with vector multiplets. When no matter is present on the brane, the background cosmology 
can be explicitly solved. In particular the metric is conformally flat implying that 

C = 0. (2.32) 

Moreover the loss parameter vanishes explicitly 

A$ 2 = 0. (2.33) 

A detailed analysis is presented in the appendix. In particular the last equation implies 
that the brane dynamics is closed. When matter density is present on the brane, we show 
in the appendix that conformal flatness is not preserved. Nevertheless we assume that the 
breaking of conformal flatness is small enough to allow one to assume that (gjgg ) and flOjj ) 
are still valid, both for the background as well as on the perturbative level. 

We now turn to the discussion of the background evolution. The solutions given below 
have been obtained in |10| |, but we need to present them in detail in order to discuss the 
evolution of density perturbations in section 4. 

2.2 Radiation Dominated Eras 

The background cosmology can be solved in four different eras: the high energy era, the 
radiation and matter dominated epochs and finally the supergravity era where the scalar 
field dynamics dominates. We discuss these four eras in turn paying particular attention to 
the matter-supergravity transition where we show that requiring coincidence now implies 
a fine-tuning of the supersymmetry breaking part of the brane tension. 

First of all we consider the high energy regime where the non-conventional p 2 dom- 
inates, i.e. for energies higher than the brane tension. In that case we can neglect the 
effective potential V and find that 



1/4 

oo ( ^ ] , (2.34) 



while the scalar field behaves like 



+ .*]„(!). (2.35) 
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In the following we will focus on the case (3 = 0, i.e. a constant scalar field, for which the 
projected Weyl tensor vanishes altogether 



£00 = 0. (2.36) 

Notice that Newton's constant does not vary in time in this case. 

The usual radiation era is not modified by the presence of the scalar field 

4> = h (2-37) 

and 

/ t \ 1/2 

a = a r {—j . (2.38) 
Newton's constant does not vary in time, while 

£00 = 0, (2.39) 
as in the high energy regime. Note that the solution in the radiation era in this brane 



world scenario is similar to the radiation era solution found in Brans-Dicke theory, see [15] 



and [16], for example. The field (j) approaches quickly the attractor for which (j> =constant. 



2.3 Matter Dominated Era 

The matter dominated era leads to a more interesting background cosmology. Let us first 
consider the pure sugra case where T = 1. This is a good approximation until coincidence 
where the potential energy of the scalar field cannot be neglected anymore. The solution 
to the evolution equations is 

<f> = <h + pln(£\ (2.40) 
a = a e (P\ (2.41) 



where r e and a e are the time and scale factors at radiation-matter equality. We are 
interested in the small a case as it leads to an accelerating universe when no matter is 
present and small time deviations for Newton's constant. 
For small a we get 

2 8 2 

7 = a . 

' 3 45 

(2.42) 

The projected Weyl tensor is given by 

4a 2 

^00 = — p-, (2.43) 
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which decreases like a 3 to leading order in a. 

In a phenomenological way we identify Newton's constant with the ratio 

8^v(r) ^ IP {%u) 
3 p m 

In terms of the red-shift z this is 

Gn(z) / Z + l\ 4a2/5 

= UnJ • (2 - 45) 

For the supergravity case with a 2 = 1/12 the exponent is 1/15. As z e ~ 10 3 this leads to 
a decrease by 37% since equality. 

Notice that the Newton constant starts to decrease only from the time of matter and 
radiation equality and is strictly constant during the radiation dominated era. Nucleosyn- 
thesis constrains the variation to be less than 20%. However, in our model we would expect 
the couplings to standard model particles to also vary in a similar manner. This could lead 
to a variation in, for example, the proton and neutron masses since these arise from Yukawa 
couplings in the standard model. We note that many of the tests for the variation of the 



Newton constant assume all other masses and couplings are constant [17|; it is possible 
that our supergravity variation would evade detection. 

2.4 Supergravity Era 

After coincidence matter does not dominate anymore; this is the supergravity era domi- 
nated by the scalar field dynamics. Let us review it briefly. Consider first the pure sugra 
case. It is easy to see that the potential vanishes 

Vsugra = (2.46) 

leading to a static universe with 

= --ln(l-a 2 £|y|) 
a 

(2.47) 

where we have defined dy = adx^. This is a flat solution corresponding to a vanishing 
cosmological constant on the brane-world. 

As soon as supergravity is broken on the brane T / 1 the static solution is not valid 
anymore. The new four dimensional potential becomes 

The time dependent background is obtained from the static solution by going to con- 
formal coordinates 

ds 2 = a 2 (u)(-d7] 2 + du 2 + dx l dxi) (2.49) 
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and performing a boost along the u axis 



a(u, r/) = a + hi], 
4>(u,r]) = (f) \ u + hrj, 



(2.50) 

where we have displayed the explicit £ dependence. Now for 

v/T 2 - 1 

^i*-^— (2.51) 

we find that the Friedmann equation is fulfilled. Similarly the Klein-Gordon equation is 
satisfied. Moreover we find that 

£oo = 0, (2.52) 

as the bulk metric is conformally flat. 

The resulting universe is characterized by the scale factor in cosmic time 

o(t) = 7T V -rj < 2 ' 53) 

with to = 3^? fer 3/2g - The scale factor corresponds to a solution of the four dimensional 
FRW equations with an acceleration parameter 

6a 2 

^rr^?- 1 <2 - 54) 

and an equation of state 

4a 2 

usugra = -1 + - - 2q2 (2.55) 

which never violates the dominant energy condition. The solution with a = — h= is 
accelerating. In particular we find that 

q = ~ (2.56) 

and for the equation of state 

5 

USUGRA = -j- (2.57) 
This is within the experimental ball-park. 

2.5 The Matter-Supergravity Transition 

Let us now investigate the transition between the matter dominated and supergravity eras. 
If we denote by Hsugra the Hubble parameter derived in the pure supergravity case, then 
the Friedmann equation in the broken supergravity case is 

H 2 = H SUGRA + "S") (2.58) 
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where we have used the fact that 4> varies slowly compared to a. The evolution coincides 
with the one obtained from unbroken supergravity as long as the contribution from the 
potential does not dominate. In the radiation dominated era this requires 



T 2k\ 3 1 - a 2 ' 

where p e is the matter density at equality. This implies that the left-hand side is much 
smaller that 10~ 39 GeV 4 . Let us now denote the supersymmetric brane tension by 



Ml = ^ (2.G0) 



3W 

and the supersymmetry breaking contribution 

M% S = {T-\)M%. (2.61) 

We find that 

M bs « (2-62) 
1 — a z 

Now this is an extreme fine-tuning of the non-super symmetric contribution to the brane 
tension. 

In the matter dominated era the supergravity Hubble parameter decreases faster than 
the potential contribution. Coincidence between the matter dominated supergravity con- 
tribution Hg UGRA and the potential energy occurs at z c 



,4 1 ( Z c +1 



3+Q/3/7 

M - = T^?lfrij (2 - 63) 

Imposing that coincidence has occurred only recently leads to 

Mbs ~ Pc (2-64) 

where p c is the critical density. This is the usual extreme fine-tuning of the cosmological 
constant. Indeed it specifies that the energy density received by the brane- world from 
the non-supersymmetric sources, e.g. radiative corrections and phase transitions, cannot 
exceed the critical energy density of the universe. 



3. Cosmological Perturbations using the Fluid Flow Approach 

We now turn to the discussion of cosmological perturbations. There are different effects 
which will influence the evolution of perturbations: 

• The evolution of the gravitational constant in the matter era changes the evolution 
of the background: In the matter era, the gravitational constant decreases and fur- 
thermore the universe is expanding slower (up to order a 2 ) than in the FRW matter 
dominated era in general relativity. 
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• Perturbations in the scalar field are the source of matter fluctuations and vice versa. 

• Perturbations in the projected Weyl tensor act as sources for the scalar and matter 
perturbations and vice versa. 

These effects will change the growth of perturbations compared to normal 4D cosmo- 
logical models or the Randall-Sundrum model. In particular it should be noted that our 
formalism allows us to treat the Randall-Sundrum cosmology. Indeed by putting a = 0, 
T = 1 and neglecting the scalar field contribution we obtain the Randall-Sundrum case 
with a flat boundary brane while putting 1 leads to a de Sitter boundary brane. 

While discussing perturbations, we use the fluid-flow approach (see |l^], [|l9| and pC|]) 
rather than the metric-based approach ||. The main difference is that all perturbation 
variables are expressed in terms of fluid-quantities, rather than metric-variables. For our 
purpose, i.e. discussing the evolution of the density contrast 5 = Sp/p of the dominant 
fluid at each epoch, this approach is simpler in order to obtain the evolution equations. To 
do so, we need to derive the Raychaudhuri equation and the Klein-Gordon equation in the 
comoving frame. 

3.1 The Raychaudhuri and Klein-Gordon Equation 

On the brane B, see figure 1, the matter energy-momentum tensor is conserved. Denoting 
by n the normal vector to the brane and defining the induced metric by 

n a b = 9ab-n a n b (3.1) 

such that n 2 = 1 and n ab n b = 0, the conservation equation reads 

V a r ab = (3.2) 

where V a = n b a D b is the brane covariant derivative and 

Tab = (P + P)u a u b - pn ab (3.3) 

is the energy- momentum tensor. Notice that the vector u a is the velocity field of the brane 
matter and thus must be orthogonal to n and satisfies u 2 = — 1. In addition we have the 
usual decomposition in terms of the shear a ab , the helicity u ab and the expansion rate 

y cU d = -u c ii d + H B u d c + a d + u; d (3.4) 

where we set a ab = 0, u) ab = later on. 
Using the relations 

V a u a = 3H B , u a = u b V b u a (3.5) 

one derives from ( |3.2[ ) 

p=-3H B (p + p). (3.6) 
Notice that no matter leaks out of the brane. Defining 

V a = u ab V b = V a + u a u b D b , (3.7) 
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which is nothing but the spatial covariant derivative on the brane, one obtains 

Ua = — , 

P + P 

whose divergence V a u a leads to the Raychaudhuri equation 

VaP 



3H B + 3H B = -R 00 - V«(— ^) 

P + P 



(3.8) 



(3.9) 



We have neglected a a b and uj a b = here. 




Figure 1: The brane- world B is perpendicular to the normal vector n. The four dimensional 
hypersurface C is orthogonal to the velocity vector u. The surface S = C n B is a set of comoving 
observers following matter on the brane. Note that the perturbed brane is not necessarily located 
at y =constant. The different metrics in the text are as follows: g a b is the full five-dimensional 
metric, h a t is the induced metric on the brane B, u a b is the induced metric on C and n a t is the 
induced metric on the comoving hypersurface S. 

In order to have a closed system of differential equations for the perturbations we 
need the Klein-Gordon equation in the comoving gauge. The metric on a hypersurface C 
orthogonal to u, see figure 1, is given by 



h. 



ab 



9ab + U a U b , 



(3.10) 
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where u 2 = — 1 and h ab u b = 0. Denning the covariant derivative on the hypersurface C by 
D a = h ab D b and using D a u a = 4H one obtains 

D a D a = D a D a - AHu c D c - u a u b D a D b . (3.11) 

In the comoving gauge we have uq = — 1 , U{ = leading to 

D 2 (p = D 2 (p-AH(p-4>. (3.12) 

We can now evaluate the Laplacian D 2 in terms of the Laplacian on the comoving surface 
S (see figure 1) orthogonal to u. Notice that V a = n ab D b . Expanding V 2 leads to 

D 2 = V 2 + (K + u a u b K ah )n c D c + n a n h D a D b (3.13) 

Now K a i,u a u b = u a D a (u • n) — u a n a , where 

K ah = D a n b (3.14) 

is the extrinsic curvature tensor. Using (u ■ n) = and n a V a p = we find that 

K ab u a u b = 0. (3.15) 

The junction conditions lead to 

K= T --U B , (3.16) 
with r = — p + 3p. Finally we can read off the Klein-Gordon equation 

4> + - 0" - (I - U B ) <P' - V 2 = (3.17) 

where prime stands for the normal derivative. Now the junction conditions lead to 

dU B 

We find that the Klein-Gordon equation reduces to 



(3.18) 



4> + 4H<P + 1^ (p - 3p) - V 2 = + A<& 2 . (3.19) 

Notice that it involves only brane derivatives when the loss parameter A<3>2 = 0. We now 
turn to the evolution equation for £"00 which enters in the evaluation of i?oo- 

3.2 The Evolution Equation for E^ v 

Let us consider the hypersurface S on the brane with induced metric u a b and orthogonal to 
the velocity vector u a such that u a i,u a = and u a u a = — 1. The induced metric is given by 
u a b = n ab + u a Ub, where n a b is the brane metric. Let V a be the brane covariant derivative 
and V a the covariant derivative with respect to u a b- 
Consider now V a E a i,: 

V a E ab = u ac V c E ab - u a u c V c E ab . (3.20) 
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The last term can be written as 

u a u c V c E ab = u c V c (u a E ab ) - u a E ab . (3.21) 
The first term in equation ( 3.20| ) can be rewritten as 

u ac V c E ab = V c (u ac E ab ) - E ab (V c u ac ) 

= V c {u ac E ab ) - E ab (u a + 3H B u a ) . (3.22) 

As a next step we consider 

V c (u ac E ab ) = V c (u ac n d b E ad ) 

= V c (u ac u d b E ad ) - V c (u ac u d u b E ad ) 

= V c (u ac u d b E ad ) - u d V c {u ac u b E ad ) - (V c u d )u ac E ad u b . (3.23) 

Using Q) 

{V c u d )u ac E ad u b = H B u ad E ad u b + (a c d + to d )u ac E ad u b . (3.24) 
Let us further define the projected tensor 

E ab = <u b d E cd (3.25) 

Then 

V c E c d = V c E c d - u c E c d . (3.26) 

So, from eq. (|3.23| ) we get 

V c (u ac E ab ) = V c E c b - u c E% - u d V c (u ac E ac u b ) - H B u ad E ad u b 

- {a d + 0J d )u ac E ad u b . (3.27) 

Collecting everything and defining E = u ab E ab ,the spatial trace, we end up with 

V a E ab = V C E% - u c E% - u d V c (u ac E ad u b ) - (u a E ab )- - H B Eu b - 3H B u a E ab 

- {a d + oj d )u ac E ad u b . (3.28) 

Let us now specialize to the comoving gauge, i.e. the surface S is a comoving one and 
u = 1, u l = 0. Notice that the first term of the last expression vanishes. We obtain 

V a E a0 = V C £ C - 4H B E 00 - E 00 - u a E% + (a d + co d )u ac E ad (3.29) 

where we have used the tracelessness of E ab to get E = Eqq . For the background we retrieve 

the previous expression as u a = and Eq = and we explicitly assume that a = 0, uj = 
for the background. 

We need further to to evaluate 

V a 7T ab = V C 7T c0 " 3# B 7T00 " H B 7t - U a TT% - 7^)0 (3.30) 

which simplifies to 

V a ir ab = V c 7r c0 - iiarf (3.31) 

and finally we need 

V a P a0 = V c P c0 - 3H B P 00 - H B P - ii a P$ - Poo (3.32) 

Collecting all these ingredients leads to the consistency equation ( [2,14 ) expressed in the 
comoving frame. 
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3.3 The Perturbed Dynamics 

We have now obtained all the necessary equations in order to derive the perturbation 
equations up to linear order. In doing so we decompose all fluid quantities into an average 
(over the comoving hypersurface) of the quantity plus a perturbation, i.e. p(x, t) = pb(t) + 
5p(x, t), and similar for the pressure and the expansion rate H. We insert these expressions 
for all quantities into the equations we obtained above and subtract the average. 

We begin with the Raychaudhuri equation and energy conservation equation evaluated 
in a comoving basis. The Raychaudhuri equation reads 

3H + 3H 2 = — — - R 00 . (3.33) 
p + p 

In this expression H is the expansion rate, p the pressure and p the energy density. All 
quantites are measured with respect to a comoving observer 1 . Rqq is the time-time com- 
ponent of the Ricci tensor. The dot represents the time-derivative with respect of the 
comoving observer (i.e. proper time). The Raychaudhuri equation expresses just the be- 
haviour of matter under the influence of geometry and is independent of the field equation. 
In particular, for matter on the brane, it is the same as in the usual 4D case. The five- 
dimensional character of the spacetime enters only through Rqq. For the FRW metric, Rqq 
is 

#00 = ~V + ^(p + 3p) + ^p (2p + 3p) - ^ 2 - E 00 . (3.34) 

The proper time is not a unique label of the comoving hypersurfaces, because it can vary 
in space. Therefore, we need to transform to coordinate time, which labels these comoving 
hypersurfaces. This transformation from proper time to coordinate time is given by (see 
§ori) 

^EL = i _ JL_ ( 3 35) 

dt p + p' v ' ' 

This gives the variation of proper time due to the perturbations. Using the conservation 

equation we obtain 

SH = H^^S - -. (3.36) 

l + uj 3(1 + a;) y ' 

where we have defined the density contrast 

5= 5 -P. (3.37) 
P 

In the rest of the paper we denote by dot the coordinate time derivative of the perturbations. 
In this paper we concentrate on the case of w = p/p = const, and c 2 s = w. We deduce that 

5H = j^-(H6 + H6) - 5 (3.38) 
1+u 3(1 + to) 

which combined with the perturbed Raychaudhuri equation 

SH = -^-H - 2H5H + - ISR 00 (3.39) 

p + p 3 a z p + p 3 



1 This observer is, of course, confined onto the brane world. 
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8 + (2 - 3lo)H6 - Quj(H 2 + H)8 = (1 + wWi^o - £<A<J (3.40) 



leads to an equation for the density constrast 

k 2 

where we have used V 2 = k 2 /a 2 . Notice that the left-hand side coincide with the usual 
four dimensional expression. The new physical ingredients all spring from 5Rqq. As the 
Friedmann equation and its time-derivative are given by complicated expressions, we do not 
derive a general equation for the evolution of 8 but rather analyse each regime separately 
and derive the corresponding equation. 

In the same fashion, we derive the perturbed Klein-Gordon equation, using eq. ( 3,35| ): 



(50)- + AH(8<t>y + 



k \pn * ,d 2 U B d 2 V 



= Us + 8$ + 4H<j>)] + [38p - Sp] . (3.41) 

1 + w L J 6 ocp 

where the coupling between the matter and scalar perturbations is explicit. 

Finally we need the perturbed version of the £"oo consistency equation. From the 
equations we derived in the previous section we obtain 

8(V a E a0 ) = V c 8E c0 - 45H B E 00 - 4H B 8E 00 - 5E 00 - -^-E m , (3.42) 

P + P 

where the shear term vanishes automatically at first order. Similarly we have 

V a Sir a0 = V c 8ir c0 - u a K« (3.43) 

and 

V a 5P a0 = V c 8P c0 - 38H B P Q0 - 3H B 5P 00 - 5H B P - H B 5P - u a Pg - (8P 00 )- (3.44) 

This allows us to write down the necessary perturbed equation. 

Notice that the components of E a Q play a role in the perturbed dynamics. As can be 
seen, they are not constrained by the brane dynamics and lead to a direct influence of the 
bulk perturbations on the brane perturbations. For this reason the perturbed dynamics is 
not closed on the brane. In the rest of this paper we will only deal with long wave-length 
phenomena where the dynamics is closed. 



4. Time Evolution of Cosmological Perturbations 

In this section we discuss some solutions to the perturbed dynamics. The equations are 
very difficult to solve in general so we restrict ourselves to the different cosmological eras. 
We start in the high-energy regime, assuming the universe to be dominated by relativistic 
particles. Then we discuss the normal (low-energy) radiation dominated epoch, followed 
by the matter dominated epoch. Because the system is not closed on the brane for large k 
(small wavelength), we present solutions only for the k — > limit. This limit enables us to 
deduce the effects of the bulk scalar field and E^ u on density perturbations. 
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In particular we can write the perturbed E^ v equation in that limit as 

5E 00 + — ^oo + ±H B 6E 00 + 45H B E 00 = -*^*p - J^E^R _ EII5 + S P 00 
p + p 4 p + p 4 4 



" P -P 00 + 3H B SP 00 + 35H B Poo + + SH B P 



P + P 



(4.1) 



Explicitly we need 



3V 30 2 - 15c/> 2 „ , 

p 00 = — + _r_ p = + _ r_ (4.2 

8 16 16 16 v 7 



The perturbed quantities are then 

SU B = U' B 5<j) (4.3) 

and 

5P 00 = y>5j> (4.4) 

as we neglect the potential in the different eras. This equation will be made more explicit 
in each era. 

As a first step we will first analyse the Randall-Sundum scenario with a flat boundary 
brane. The differential equation for the density contrast is of third order leading to the 
appearance of three modes. This is to be compared to the usual two FRW cosmology modes 
derived from four dimensional general relativity. As already mentioned two of the modes 
will coincide with FRW cosmology while the third mode is entirely due to perturbations 
in the Weyl tensor signalling a breaking of conformal invariance in the bulk. 

4.1 The Randall Sundrum scenario 
4.1.1 The high energy era 

In the high energy regime, the equations for 5Eqo and 5 read 

5 + H5- im 2 5 = -Ue o, (4.5) 

{8Evoy+4H{8Ew) = Q. (4.6) 

This leads to 

SE 00 = (4.7) 

The overall solution of the equation for the density contrast is a sum of the solution of the 
homogeneous equation as well as the general solution. It can easily be found as 

6 = 6o& 2 + 6 1 t- 3 '* + ^6E t. (4.8) 
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4.1.2 The low energy radiation era 

In this era, the equations for ^-Eoo and $ read 

5 + H5-2H 2 5 = —SEm, (4.9) 
(5E 00 )- +4H(6E 00 ) = 0. (4.10) 

Again, the solution to the last equation is 

ftE 

5E m = (4.11) 
The full solution to the first equation is found to be 

5 = S t + S 1 t- 1 / 2 + ^SE . (4.12) 

We do, therefore, find the normal growing and decaying modes in this regime and a constant 
mode, which is absent in FRW cosmology. 

4.1.3 The matter dominated epoch 

Finally, in the matter dominated epoch the equation for 5Eqq and 5 read 

6 + 2H 6 - ^H 2 5 = -5E 00 , (4.13) 
{5E m y + AH{5E m ) = 0. (4.14) 

The solution to the last equation is 

5E m = 5E t- 8 / 3 , (4.15) 

and the solution to the first equation can then found to be 

5 = 8t 2 l* + 5 l t- 1 + 9 -5E t- 2 l\ (4.16) 

We recover thus the normal growing and decaying mode as well as a mode sourced by 5Eoo, 
which is decaying rapidly. 

We now turn to the case with a scalar field in the bulk. The previous Randall-Sundrum 
modes will be modified in two ways. First of all the brane potential will lead to slight 
deviations of the mode exponents similar to the modification of the scale factor exponent 
in the matter era. Then there will also appear new modes due to the fluctuations of 
the scalar field governed by the Klein-Gordon equation. We will now analyse each of the 
different regimes in turn. 

4.2 Effects of the bulk scalar field on cosmological perturbations 
4.2.1 The high energy regime 

In this regime, the terms which are quadratic in p and p dominate both in the background 
and in the perturbation equations. We assume, that relativistic particles dominate the 
expansion, i.e. p = p/3. 
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Considering the large wavelength limit is sufficient because all cosmologically relevant 
scales are far outside the horizon. The perturbation equations in this regime are 

5 + HS = 18H 2 5- -5E 00 , (4.17) 
(S<j>)" + 4H(6<l>y = 0, (4.18) 
(SE QO y + 4HSE Q0 = -^M{S<j,y (4.19) 

where we have assumed that 3H(5<p)' 3> VScf). 

To obtain the solutions to these equations, we first consider the Klein-Gordon equation. 
The solutions are 

5<f> = Sfo + Sfahit, (4.20) 

where 5(f>o is a constant mode. We can now find 5Eqq. 

6E 3aU B ,, hit 
SE 00 = — —Hi— (4.21) 

Notice that the solution comprises two parts. The \jt mode is a solution of the homoge- 
neous equation while the logarithmic mode solves the complete equation. In the following 
we will always find that the solutions are expressible as a sum of homogeneous modes and 
modes solving the complete differential equations. 
We can now deduce the density constrast 

5 = 5 t 3/2 + 5!t~ 3 / 4 + 5 2 (t) (4.22) 

where the complete solution reads 

S 2 (t) = -U*/ 2 j l dt'(t')- 13/4 f dt"5E 00 (t") 7 / 4 (4.23) 
In the long time regime, we focus on the leading growing mode, obtained by approximating 

5 2 {t) = 0((lnt)t) (4.24) 

This implies that we find two leading growing modes in i 3 / 2 and (lnt)t. Notice that the 
logarithmic mode is triggered by the scalar fluctuation, <50i, and is therefore absent in 
the Randall-Sundrum case. Moreover for very long times we see that the leading mode 
increases like i 3 / 2 , which is much larger than in FRW cosmology. Interestingly, this growth 
of fluctuations is anomalous, in the sense that the exponent of the growing mode is larger 
than one. The reason for this is the source term (18H 2 5) which is much larger than the 
normal 4TrGp m 8. In addition, there is the contribution of bulk gravity. In this regime we 
cannot expect the normal behaviour, because gravity is simply not four-dimensional. 

4.2.2 Radiation domination 

In this regime the Newton's constant is not varying in time. Moreover matter does not 
appear in the Klein-Gordon equation. This leads to the perturbation equations for the 
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scalar field and the density contrast 

(6(/>)"+4H(6<t>y =0, (4.25) 

5 + H5-2H 2 5 = 4aH 2 5<p-UEoo, (4.26) 

(6E 00 )- +4HSE 00 = -3aH 2 (5(t)y. (4.27) 
As before we find homogeneous and complete solutions. More specifically 

6<f> = 5(j) + 5(f)-t- 1 , (4.28) 

leading to the perturbed 5Eqq 

. ,_, 5Eq 3a . , In t . A 

5E 00 = ^ + —50--p. (4.29) 

Notice that the first mode springs from the homogeneous equation and the last one from 
the complete equation. The density contrast can then be deduced: 

S = Sl t + <5_ 1/2 rV2 + t j l dt'(tT 5/2 f (t^W ~ |(f ) S/a «*») (4.30) 

There are two homogeneous modes as in FRW cosmology. New contributions emerge from 
the scalar field. In particular there is a growing mode in 0(54>- kit) which is triggered by 
the decreasing scalar mode in 8(j>-/t. 

4.2.3 Matter domination 

In this regime Newton's constant varies in time as the background scalar field is time 
dependent. This leads to a very rich structure of modes for the density contrast. The 
perturbation equations read 



6 + 2H5 = ^aH 2 



a 



SE 00 , (4.31) 
(6<f>)" + AH{5^)- + 2a 2 H 2 {5(f>) = -2aH 2 5, (4.32) 



(5E 00 y + 4H(5E 00 ) + 4(5H)E 00 = -3aH 2 (5<py - 3aH 2 j5 

+ ^ [-t" 2 W + t- 1 W] + 3H^(S<t>y - \pH-H. (4.33) 

This system of equations possesses power law solutions, which we derive up to order a 2 . 
Defining each mode by the ansatz 

6 = Sit a > , 5<p = 5<pit b * , 5E 00 = 5E l 00 t c * (4.34) 

we find that 

at = bi = <h + 2 (4.35) 

for all modes. There are two types of modes. Let us first discuss the complete solutions 
of the coupled differential equations. We find that there are three modes corresponding to 
the exponents 

2 11624 o 1036 2 2 808 2 /A , 

^ = 3 " W° ' = - 1 " ^75 ° ' a ~ 2 '* = -3 + 225 a (4 " 36) 
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The modes #2/3 and S-i are deformations of the FRW modes due to the scalar field. In 
particular we find that a 2 /3 ~ 0.55. Thus, the growth of fluctuations is smaller, than in the 
normal matter dominated epoch. The perturbation of the bulk scalar field is growing with 
the same exponent as the density constrast and SEqq is rapidly decaying with an exponent 
c = a - 2 Pt -1.45. 

In addition to the growing mode, there are two decaying modes. One of these decaying 
modes is a modification of the usual a = — 1 decaying mode. We obtain a = —1.23. Thus, 
the mode is decaying even faster than normal. The other decaying mode corresponds to 
a = -0.37. 

On top of these three modes there are seven homogeneous modes corresponding to 
zero modes of the various differential operators appearing in the three coupled equations. 
Let us discuss them in some detail as they have remarkable properties. The homogeneous 
modes of the perturbed Raychaudhuri equation have characteristic exponents 

fl2 / 3 = 3 + 225 a ' a - 1 = " 1 + 225 a (437) 

Notice that these modes are two deformations of the usual FRW modes. The scalar field 
potential lifts the degeneracy between the doublets of a = 2/3 and a = — 1 modes. We find 
that a 2 /3 ~ 0.67 very close to the FRW value. That the growing mode is larger than the 
one found in general relativity comes from the fact that the matter dominated universe 
is expanding slower in the brane world model. Therefore, overdensities, which want to 
contract but have to compete with the cosmological expansion, can grow more easily in 
the brane world model. 

The homogeneous modes of the Klein-Gordon equation are characterized by the expo- 
nents 

82 ^ ^6 2 , . 

ao = -—a , a_ 5/3 = -- + —a , (4.38) 

corresponding to an almost constant mode and a mode which decreased very fast. 

The E 00 equations comprises three differential operators. The first order differential 
operator acting on 8Eqq leads to a single mode with 

2 32 

a-2/3 = - 3 + ^« 2 - (4-39) 

The differential operator acting on 8(f) leads to two modes 

512 

a = 0, a_i = -10 + —a 2 . (4.40) 
Eventually the operator acting on 5 leads to 

10 

«o = ^. (4.41) 

Notice that in the epoch of accelerated expansion the structures get frozen in. We have 
thus described the perturbation modes in the four cosmological eras. 

In conclusion, we have found two the growing modes, with exponent a 2 /3 ~ 0.67 
and a 2 /3 ~ 0.55. The 00-component of the Weyl-tensor is always decaying, whereas 
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perturbations in <f> follow the density contrast. It is instructive to compare these results 
to the ones found in Brans-Dicke theory. Here the exponent of the growing and decaying 
modes are modified as follows pil: 



2 2 1 

a+ = ~ + — and a_ = -1 - — , (4.42) 

3 3ll> 3lo 

where oj is the Brans-Dicke parameter. Note, that there are more modes in these theories 
as well. Here we find that the growing and decaying modes are shifted due to the presence 
of a scalar field which couples to gravity in a way which differs from Brans-Dicke theory. 

4.2.4 Observational consequences 

Although we have focused on the large wavelength limit only, we can draw some (quali- 
tative) conclusions from our findings. First of all, which of the two growing modes in the 
matter dominated epoch appear, depends on the initial conditions imposed in the early 
universe. It is likely that both modes are generated. If the mode with exponent a ~ 0.67 
is not generated and the power spectrum is normalized to the observations today (such 
as to galaxy clusters), the brane world model would have much more power in the matter 
perturbations than one would expect in normal quintessence models, for example. Thus, 
in this case there should be more galaxy clusters at high redshift in the brane world model 
than in normal Einstein gravity with comparable matter density (with a cosmological con- 
stant or quintessence field). If the mode with exponent a ~ 0.67 contributes significantly 
to perturbations, then the power in the perturbations in both theories would be similar. 
However, much more work is needed in order to make this more concrete. 

Another consequence is a modification of the spectrum of anisotropics in the cosmic 
microwave background radiation (CMB). This is because the distance to the last scattering 
surface (LSS) will be modified due to the slower expansion in the matter dominated epoch. 
As a result, the first peak will be shifted to larger angular scales. However, we believe 
this effect to be small to be detected. In addition to this effect, it is conceivable that 
isocurvature modes (between the Weyl fluid and radiation, for example) might survive 
quite long in the radiation era, leaving their imprint in the CMB. To make more concrete 
predictions, it is necessary to go beyond the small fc-limit, which we considered in this 
paper. This involves a study of the bulk gravitational field, as well as perturbations in the 
bulk scalar field away from the brane |Q] . The calculations are very difficult and a detailed 
discussion of the bulk equations is beyond the scope of this paper and will be presented 
elsewhere. 



5. Conclusions 

We have presented a discussion of cosmological perturbations in brane world scenarios with 
a bulk scalar field. For a model motivated from SUGRA in singular spaces, we were able 
to find solutions of the perturbation equations in the large scale limit. As is the case with 
brane world scenarios of the Randall-Sundrum type, the pertrubation equations on the 
brane are not closed. Instead one has also to solve the bulk equations, which is a very 
difficult task in general. 
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Nevertheless, our findings indicate that there are considerable differences to the Randall- 
Sundrum scenario as well as to usual scalar-tensor theories in four dimensions. The growth 
of structures is different than in these models. A more elaborate investigation of perturba- 
tions, in particular of the bulk perturbations, would be necessary in order to make definitive 
predictions concerning cluster abundance, large scale structure and the anisotropies in the 
CMB. 



There are other open questions, which our work leaves. As discussed in [22], in models 
based on SUGRA in singular spaces we expect a second brane in the bulk at some distance 
from our brane universe. We have not addressed the influence of this mirror brane on 
the structures on our brane universe. Similarly, we haven't discussed the dynamics of the 
radion in these models. The effects of these should be encompassed in the Weyl tensor. 
We will, however, turn to these questions in future work. 
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A. Supersymmetric Backgrounds 
A.l BPS configurations 

In this appendix we will consider configurations which preserve supersymmetry in the bulk. 
To do so we will look for Killing spinors satisfying the identities 

S e ^i = 0, 5 e Xf = (A.l) 

where is the gravitino spinor (i = 1, 2) and A? belongs to the vector multiplet comprising 
the scalar field <p x . This leads to the first order equations 

D^i + % -l»WQ l3 e j = (A.2) 

o 

and 

i-fd^f* + \Y r (hj< ! = (A.3) 
where Q\ = Qa{& a )l and Q a Q a = 1- Define the vector p a such that 

and the matrix 

r 

such that r 2 = 1. Rotation invariance implies that the only non-zero components are 
p z = 1 and pt- The unknown component pt is fixed by the boundary condition. Now the 
spinors can be split into positive and negative chiralities 



--Pad z 

Pal" 1 



(A.4) 
(A.5) 



(A.6) 



satisfying 

ef = : iYQ;,< ■ ' (A.7) 

At each point this select the chirality defined by T, i.e. we preserve one half of the super- 
symmetries. Now ( |A.3[ ) is satisfied provided 

d z <t> x = ±^=W' X (A.8) 



This is the BPS equation for supersymmetric configurations. Now choosing the positive 
sign, we find that the boundary condition at z = is automatically satisfied provided 

Tv 7 ^ = 1 (A.9) 

This leads to the vanishing 

A$ 2 = (A.10) 
Let us come to the bulk evolution of the Killing spinors. They satisfy 

D^i + '-l^WQ^ = (A.ll) 

o 

This is a first order differential equation which can be solved provided the integrability 
condition 

[D a , D b ]ef = Rabcd^ef (A. 12) 



is fulfilled. Now one gets 



[D a ,D b ]et = -^WSaO + l Palb ™ 7a W' x W, x T)et (A.13) 



from which we deduce that 

Kbcd = -j^W 2 (g ac g bd - g a d9bc) + -^W^ x W iX (p a g b[d p c] - p b g a [dPc}) ( A - 14 ) 

For W )X = 0, i.e. at the critical points of the superpotential one recognizes the Riemann 
tensor of AdS^. More generally we find that the background geometry which preserves 
supersymmetry in the bulk is such that the Weyl tensor vanishes 

W abcd = (A. 15) 

as the Weyl tensors form, in the set of curvature tensors, the complement to the antisym- 
metrized product of the metric tensor with itself or with a symmetric tensor. This implies 
that the bulk geometry is conformally flat and can be written as 

ds 2 = e 2A ^ t] {-dt 2 + dz 2 + dxidx 1 ) (A.16) 

This leads to the vanishing of the projected Weyl tensor 

E ab = (A.17) 
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in the bulk. Moreover the bulk and brane expansion rates coincide 



H B = H. (A.18) 
Notice that ( A.14| ) leads to the Ricci tensor 

W 2 W' X W X 2 

Rab = — Y 9ab H 42' ^ 9ab + 3 ^Pfe) (A.19) 

and the curvature scalar 

hW 2 

a 2 R = _<^_ + 2W x Wx (A.20) 
This implies that Einstein equations 

Rab ~ ^9ab = T ab (A.21) 

are satisfied with 

d z A = — ^ (A.22) 

Notice that ( [A. 22 ) and (A.8) are the BPS equations which lead to the accelerating universe. 
Solutions of these equations satisfy the Einstein equations and the Klein-Gordon equation. 

Let us now consider the second brane where supersymmetry is not broken. Now the 
Killing spinors for the super symmetric T = 1 case satisfy the projection equation 



ef = ±h 5 Qije^ (A.23) 



as pt = 0. This is not compatible with ( A. 7 ) as T 7^ 7 s . The only solution is therefore 
e = 0, i.e. supersymmetry is completely broken by the non-supersymmetric brane with 
T / 1. Notice that the breaking of supersymmetry is global due to the presence of two 
boundaries respecting incompatible supersymmetries. 

So we have shown that breaking supersymmetry on the brane leads to broken = 
background configurations which still satisfy a system of two first order BPS conditions. 

A. 2 Breaking Conformal Flatness 

Let us now consider the case where matter and radiation are present on the brane. We will 
show that one cannot deform the bulk geometry in such a way that the boundary conditions 
are satisfied. This implies that matter on the brane breaks the conformal flatness of the 
bulk. 

Let us perform a small change of coordinates in the bulk 

x = x - £ (A.24) 
inducing a variation of the bulk metric is 

S9ab = d a Cb + d b i a (A.25) 

and of the scalar field 

5(3 a 0) = d^da?, 5{d a d b <t>) = d c 8 a 0d b e + 3 c c> b <R£ c (A.26) 
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We also need to evaluate the variation 



S(d z g ab ) = d z £ c d c g ah + d a £ c d z g cb + d b ^ c d z g ca (A.27) 
The boundary equation 5(d z (j))\o = leads to 

d z e\o+Ptd z e = (A.28) 

We then find that 

S(d z9aa )\o = 2d a ^d z g aa \ (A.29) 

with no summation involved. The metric boundary condition at the origin are modified 
according to 

r / ''z 9aa \ i aa) I @z9aa i aa i [ a or .\ 

<H )|o = |o lo lo (A.30) 

9aa 9aa 9aa 9aa 

implying that 

5(^)\ = -^d a g aa £ a \o (A.31) 

9aa 2 

with no summation involved. The ga boundary condition then reads 

^ Vfclo = ~ (A.32) 

As the background is x l independent this cannot be satisfied unless p = 0. This proves 
that the bulk metric cannot be smootly obtained from the matterless case by performing 
a change of coordinates in the bulk. In particular conformal flatness is broken leading to 
an explicit breaking of supergravity by matter on the brane. 
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